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■ Abstract 

■ Hadronic decays of B mesons are reviewed. First, masses of B mesons and ob- 
served patterns together with physics behind them are discussed. Then the effective 
Hamiltonian responsible for major decays is presented and its practical applications 
are discussed in the context of factorization. Various tests of factorization are 
then studied. For rare decays, the focus is placed on Kit, tttt final state and the 
penguin-mediated X s ^. In general, the measurements are in excellent agreement 
with predictions of the standard model. 

(N 

■ 1 Basic Methods on Upsilon-4S Resonance 

On ■ 

Most of the data presented in the following are collected on the upsilon-4S resonance, and 
1 \ some basic experimental techniques are briefly described below. 

The B meson pair production cross section on the upsilon-4S resonance is roughly 
1 nb; namely, an integrated luminosity of 1 fb _1 would generate 1 million B meson pairs. 
The CLEO-II detector has logged about 1.2 fb _1 of data thus generating 1.2 million B 
^ ■ meson pairs. 

c3 ' On the upsilon-4S resonance, light quark pairs (uu, dd, ss, and cc - often referred 

to as the 'continuum') are also generated in addition to the B meson pairs. The cross 
section ratio of B meson pair to the continuum is roughly 1 to 2.5. The continuum is 
often a major background and in order to understand this component, data are taken 
right below the resonance (32 MeV below the peak) corresponding to about one half of 
the integrated luminosity taken on the resonance. When we want to plot a distribution 
of certain parameter for B meson pairs, we can subtract the distribution for the data 
taken off-resonance from that taken on- resonance (with a proper normalization). The 
distribution is then said to be 'continuum subtracted'. 

At the upsilon-4S resonance, the B mesons are generated with definite energy and 
momentum given by 

E B = £ bC am = 5.289GeV, P B = 0.325GeV/c (1) 



*A talk given at 'B Physics at Hadron Colliders', Snowmass 1993. 
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When reconstructing a decay B — > f\ + / 2 H — • / n , natural parameters to look at are thus 
the total energy and momentum of the decay products (i = 1, ..,n): 

E tat = Y,Ei, Ptot = Y.Pi ( 2 ) 

which should peak at Pbeam and Pb respectively, where E± and Pj are the energy and 
momentum of the i-th decay product. In practice, often used parameters are the 'energy 
difference' AE and the 'beam-constrained mass' Mb defined by 

AE = E tot - P bcam , M B = vX 2 cam - Pit- (3) 

Since Pbeam is a constant, measuring AE and Mb is equivalent to measuring E tot and 
P tot . The mass reconstructed this way has a good resolution which varies from 2.5 to 3.3 
MeV depending on decay mode and usually dominated by the spread of beam energy. 
The essence of this method in background rejection, however, is simply the conservation 
of energy and absolute momentum in a B meson decay. We will often be referring to M B 
and AE in the rest of this article; the definitions are as defined above. 

2 Masses 

2.1 B andB° 

The masses of neutral and charged B mesons can be measured by fully reconstructing the 

major decay modes. Figure ffl shows the distribution of the beam-constrained mass Mb for 
— o 

B~ and B mesons after requiring that the energy difference AE is within 2.5<r of zero. 
The decay modes used are B~ — > D*°tt~ , D*°p—, D°tt~ , D°p~ ,ipK~ for the charged B 
meson and B° — > D* + n~, D* + p~, D + n~, D + p~,ipK*° for the neutral B meson. The D* 
mesons are detected by the decays D*° — > D°7r°, D* + — > D°tt + and D mesons are detected 
by D° — > K~n + ,D + — > K~7t + tt + . These modes are chosen since they are particularly 
clean. There are 362 signal events for B~ and 340 signal events for B . With a correction 
due to initial state radiation of -1.1 ±0.5 MeV, we obtain M B o = 5280.3 ±0.2 ±2.0 MeV 
and M B - = 5279.9±0.2±2.0 MeV . The first error is statistical and the second systematic. 
The systematic error is dominated by the uncertainty in the energy scale of the storage 
ring which cancels when we take the mass difference: Mbo—Mb- = 0.44±0.25±0.19 MeV; 
namely, the masses of B~ and B are consistent with being identical within several tenth 
of MeV. The results are summarized in Table [I] together with previous measurements. 

It is interesting to compare this result with that for strange and charm mesons. 
There we have M K o - M K - = 4, 024 ± 0.032 MeV and M D+ - M D o = 4.77 ± 0.27 MeV || 
which seem to indicate that the meson mass is heavier when a heavy quark is combined 
with a d quark than with a u quark. The pattern, however, clearly does not repeat for B 
mesons. The current understanding for the isospin mass splitting is that there are effects 
due to the u — d mass difference as well as QED effects || (i.e. due to the electric charge 



2 



(MeV) 


CLEO 1.5 


ARGUS || 


CLEO II § 


M B o 
M B ~ 
M B o - M B - 


5278.0 ±0.4 ±2.0 
5278.3 ±0.4 ±2.0 
-0.4 ±0.6 ±0.5 


5279.6 ±0.7 ±2.0 
5280.5 ± 1.0 ±2.0 
-0.9 ±1.2 ±0.5 


5280.3 ±0.2 ±2.0 
5279.9 ±0.2 ±2.0 
0.44 ±0.25 ±0.19 



Table 1: Masses of neutral and charged B mesons. 



Figure 1: The beam-constrained mass for charged (a) and neutral (b) B mesons after AE 
is required to be consistent with zero. Particularly clean modes are selected and summed. 

difference between u and d quarks). Both are of order a few MeV, and the two kinds 
of effects happen to cancel for the B meson case |J. There seems to be no simple and 
intrinsic reason to give M B = M B . 

2.2 Other Bottom Mesons 

Bottom hadrons heavier than B~~ and B~° are not produced on upsilon-4S resonance, and 
the results so far come from accelerators that operate at higher energies. 

Figure ^| shows the decay B$ — > ip(p, <fi — > K + K~ observed by the CDF collaboration 
in pp collisions at 1.8 TeV cm. energy. There are 14 ±4.7 events observed and fitting 
a gaussian to the peak, the B$ mass is determined to be 5383.3 ± 4.5 ± 5.0 MeV. The 
ALEPH collaboration has also reported a result on B$ mass from two events B$ — > ip'4> 
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Figure 2: B s — > ^0 decay observed by the CDF collaboration, (a) Invariant mass of 
il)K + K~ when the K + K~ pair forms the mass (within ±10 MeV). The dots are for a 
mass side band, (b) Invariant mass of K + K~ when the <pK + K~ mass is in the B s peak 
(within±20 MeV). 
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Modes 


Number of events 


M Bs (MeV) 


ALEPH § 
CDF § 
OPAL0 

DELPHIQ 


Dlj(iT-orai),ip(j) 


2 

14 ±4.7 
(1 candidate) 
4 


5368.6 ± 5.6 ± 1.5 
5383 ±4.5 ±5.0 
5360 ± 70 
5357 ± 12 ±6 



Table 2: Measurements of Bg meson mass. The ifiO and <fi mesons are detected by 
ip( ) — > l + l~ and — > K + K~ respectively, and mesons are detected in the modes 
Dg — > 0vr + , K* K. 



and Dg7r~ . The mass measurement is dominated by the event and gives 5369±5.6±1.5 
MeV. These results are summarized in Table |^ together with a possible candidate event 
reported earlier by the OPAL collaboration and a recently reported result from DELPHI. 
The measurements by CDF and ALEPH are marginally consistent (2-sigma difference 
statistically); taking the weighted average, the mass difference between Bg and B° is 97 
MeV. The value is strikingly similar to the charm case M D + — Mp+ = 99.5 ± 0.6 MeV [I], 
and also consistent with predictions of non-relativistic models: M Bs = 5345 — 5388 MeV 
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The mass of B*(J P = 1") has been measured by CUSP @ and CLEO @ by 



detecting the monochromatic photon in the transition B* — > B"f. The numbers are |T4] 

Mb* -M b = 46.4 ± 0.3 ± 0.8 MeV (CLEO) (4) 
45.6 ±1.0 MeV (CUSP). (5) 

These measurements are in accordance with an intriguing observation on the hyperfine 
splitting 

AM = M 2 (r) - M 2 (0") = const w 0.5GeV 2 . (6) 

This holds well for (vr,p), K, D, Dg and now for B. In non-relativistic models, such 
relation is realized when the potential between the constituent quarks is linearly increasing 
as a function of the distance between the quarks [[□], p^5| . It is consistent with a naive 



picture that the two constituent quarks are connected by a flux tube with a constant 
tension. At short distance, the potential is expected to be Coulomb-like; this portion 
of the potential, however, is not expected to play a significant role ]T3. Also, there is 



an electromagnetic hyperfine splitting which violates the relation |6], but its effect is also 



much smaller than the hyperfine splitting due to strong interaction [17 . 

Apart from the theoretical importance, the above mass difference indicates that B* 
cannot decay to Bit. It has a practical implication that one cannot tag the sign of the 
bottom flavor by the decays such as B* + — > B°tt + where the charge sign of the pion 
tells us if the neutral B meson is bottom or anti-bottom. Such flavor tagging would have 
made it easy to study the CP violating decay asymmetry in B° or B — > ipKg, 7r + 7r~ 
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etc. particularly in hadron colliders. Now we have to hope that there may be a higher 
resonance that decays to Bit which is narrow and produced copiously [p"8|l . 



3 Non-Suppressed Decays 

3.1 Effective Hamiltonian 

The interaction of interest for B meson decays comes from the charged current part of 



the Standard Model Lagrangian [19 



V2" 



d 
s 
b 



C C c = -^={u,c,t) Lltl V\ s I W. (7) 



L 



where g is the weak coupling constant, the subscript L for the quark field indicates 
left-handed component (e.g. Ul = |(1 — lb) u etc.), and the matrix V is the Cabbibo- 
Kobayashi-Masukawa (CKM) matrix: 



V 



' Vud Vus Vub 

V cd Vcs V cb I . (8) 

V v td v ts v t 



lb 



The experimental value of the CKM matrix V is well represented by [^, ^TJ (assuming 
unitarity of V) 

A ~ 0.22 

V ~ I -A 1 A 2 j where { a = arg(Kb) (9) 




(3 = aig(Vtd) 

and the magnitude of V u b,Vu is of order A 3 . Taking the first and third columns, the 
unitarity condition 

Kd v ub + v c * d v cb + v; d v tb = o (io) 

becomes a triangle as below (called the unitary triangle). 

(11) 
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At energy scales well below the W mass, the propagation of W can be 'integrated 
out' and we obtain 4-fermion effective Hamiltonian [22] relevant to B decays given by 



:i2) 



Q 1 = (du)(cb), 2 = (cu)(db) (13) 

where Gf = g 2 / (i\/2M^) is the Fermi coupling constant and the quark current (q'q) is a 
short hand for g^7 M (l — '~f5)q a which is a color-singlet V — A current (a: color index). Any 
combination of replacements c — > u, u — > c and d — > s can be made to obtain other possible 
interactions as long as the replacements are consistently made including the indexes of 
the CKM matrix elements. The terms shown in ([12]) are part of an expansion of the 
effective hamiltonian (the operator product expansion p3| ). It has an advantage that 
the calculable short-distance effects are separated into the coefficients of the operators 
(Wilson coefficients) while the long distance effects such as the state of valence quarks in 
mesons are absorbed into matrix elements of the operators. 

Without QCD correction, we only have the first operator Oi which is shown dia- 
grammatically in Figure 0(a). With QCD correction, gluons flying between the quark lines 
can shuffle the color flows and generate an effective neutral current operator O2 shown in 
Figure |(b). The Wilson coefficients can be calculated using the leading- logarithm 
approximation (LLA) 



Ci 



-(C+ + eg 



c 9 



CL 



with 



where oL 
given by 



■2 b 



(14) 



(15) 



—2d + = 8, and as is the running coupling constant of strong interaction 



An 



with b 



11 - -n f . 

3 ; 



(16) 



Mog(/i 2 /A^ CD ) 

with rif being the number of relevant flavors, and /1 the typical mass scale of problems 
in question. Note that C + and C_ are related by C\C- = 1. With \i 
rif = 4, and Aq CD = 0.25 GeV we have 



5 GeV, 



Ci(m 6 ) = 1.11 C 2 {m b ) = -0.26. 



(17) 



The next-to-leading logarithm approximation (NLLA) has been computed ; the result 
does not differ drastically from the LLA result quoted above. For the transition b — ► esc, 
however, the momentum transfer associated with the light quarks are much smaller than 
the bottom mass scale and as a result the corresponding coefficients could be significantly 
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Figure 3: Four fermion operators of the effective Hamiltonian responsible for B meson 
decays. 



different from flT^). In fact, in one estimation using heavy quark effective theory (HQET) 



| 3lfl , the coefficients are about 30% larger for C\ and almost twice as large for C 2 j32| : 

Ci ~ 1.45 C 2 0.45 (for b -> esc) (18) 

There are also 4-fermion operators of the type shown in Figure 0(c) called Penguin 
operators | 26fl . The corresponding coefficients, however, are small and the Penguin oper- 



ators are relevant only for highly suppressed decays such as B — > K*^y and Kir, to which 
we will come back later. 

3.2 Two-body Decays and Factorization 

Compared to semileptonic decays, hadronic decays are harder to understand due to variety 
of short and long-distance strong interactions among the quarks involved. Two-body 
hadronic decays, however, are the simplest kind, and some framework of understanding 



- factorization - exists p7| . Also, it should be noted that two-body decays account for a 
substantial fraction of total hadronic decays of heavy mesons (~ 15% for bottom mesons 
and ~ 75% for charm mesons when resonances are included |28fl ). 

The idea of factorization for hadronic weak decay dates back at least to the early 60's 
when Schwinger showed that the AI = 3/2 transition of K — > tctt can be estimated from 
the corresponding semileptonic rate [p9| . The procedure, however, was not considered to 
be accurate; in fact, when Feynman reported calculations of A — ► pir and K + — > 7r + 7r° 
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Figure 4: Decay B — > D + 7r by the operator Oi (a) and 2 (b). The latter is suppressed 
by a factor £. 



using the idea of factorization |30 , he preceded the discussion by the following disclaimer: 



'You may not wish to consider this line of flimsy reasoning; we are becoming very uncertain 
about this matter, nevertheless I shall present it.' There is, however, a good reason to 
believe that the factorization works well for certain B decays. 

We take B — > D + n~ as an example. This can occur by the operator Oi as shown 
in Figure ^ (a), where it is assumed that the B — ► D transition is caused by the current 
operator (cb) and that tc" is created by the current operator (du). Assuming that the 
B — > D transition and the 7r~ creation are independent, the amplitude can be written as 

(D + n-\(du)(cb)\B°} = (7r-\(du)\0}(D + \(cb)\B ) (19) 

which constitutes the essence of the factorization assumption. 

It is instructive to visualize the situation intuitively. A B meson may be viewed as 
an analog of a hydrogen atom where the heavy bottom quark is at the center surrounded 
by a cloud made of light quark and gluon [Figure 0(a)]. Upon the decay of the b quark, 
the b quark disappears and c, u, and d quarks appear. The c quark will combine with 
the original cloud that was around the b quark to form a D meson, and the ud pair will 
eventually turn into a pion. Here one can cast doubts on the factorization assumption on 
two points: 

1. When the ud pair passes through the cloud, it may strongly interact with the cloud, 
in which case the formation of the D meson and the creation of the pion cannot be 
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Figure 5: An intuitive picture of the decay B — > D + ir . Before the decay (a), imme- 
diately after the b quark decay (b), and right after the formation of final state mesons 
(c). 



independent. 

2. After the D meson and the pion are formed, they may re-scatter through final-state 
interaction (FSI); e.g. D + + 7r~ — > D° + 7r° etc. 

For each of the above, Bjorken has argued that it does not pose serious problem for the 
factorization assumption f33| . First, the invariant mass of the ud pair is of order pion 



mass; thus, they are highly collinear and close together. Since the total color of the pair 
is zero, they form a small color dipole and the cloud cannot see them from some distance 
away. The pair is thus expected to pass through the cloud without much interaction. 
Second, the formation time of the pion in its own rest frame is of order 0.3 fm/c which is 
the time for light to propagate from the center of the pion to the edge. Since the pion is 
highly energetic (~ 2.5 GeV), by the time it is formed the distance between the D meson 
and the pion is already several fermis; thus, they cannot interact through FSI. A similar 
argument of 'color transparency' was also used for production of p and ip in high energy 



scatterings [34 



This line of argument has been put forward by Dugan and Grinstein in the framework 
of QCD and the heavy quark effective theory, and it has been shown that factorization 
holds in the limit of Mb,d oo while Mb/Md is kept constant |35| . For decays which 



involve two charmed mesons such as B° — ► D S D + , the two mesons in the final state are 



10 



partially overlapped at the formation time, and thus the factorization may not work well 
for these decays. Factorization is known to hold also for the large Nq limit where Nq is 
the number of colors H36| . Even though the correction to the limit is of order 1/3 which is 
quite large, the applicability of the 1/Nc argument is not restricted to the large velocity 
limit ||37|| , and thus complementary to the 'color transparency' argument. 

The decay B — > D + n~ can also proceed by the operator O2 as shown in Figure £|(b). 
In this case, naively only the color singlet component of the u and d legs is expected to 
contribute. Applying Fierz transformations to color indexes as well as to gamma matrices 



|38fl, O2 can be written as 

2 =^0 1 + ^(d\ i u)(c\ i b) (20) 

where the second term is a color singlet operator formed by two color-octet currents 
with A* being the SU(3) Gell-Mann matrices. Thus, O2 contains Oi within itself, and 
consequently Oi and O2 are not orthogonal [57]. The overall coefficient of Oi is then 



C\ + C2/3. For the decay B° — > D°ir°, the relevant operator is O2. There, the role of Oi 
and O2 are inverted with the overall coefficient of O2 being C 2 + C\/3. In fact, we can 
write (0) in two ways 

C 1 1 + C 2 2 ={C 1 + f)0 1 + \(dXu){c\ l b) 

= (C 2 + f )0 2 + \(d\%)(c\ lU ) 1 > 

Assuming factorization, the effective Hamiltonian may then be written in terms of 'fac- 



torized hadron operators' |39| as 



Gf - - 

Whad = -^KdK&M rfM )had(c&)had + a 2 (db) had (cu) haid ] (22) 

where the above arguments suggest 

ai= Ci + £C 2 1 ( . 

a 2 = Ci+tCi Wlth e= 3' (23) 

where the effect of O2 to the first term and that of Oi to the second term is parametrized 
by £ (sometimes called 'color suppression factor'). The contribution of the octet current 
term in (|2~0D , however, may have a significant effect; in fact, an estimation using QCD 



sum rule indicates that its contribution may in effect lead to £ ~ El . Also, an analysis 



of charm decays suggests £ near zero p8| . It has thus been suggested that a%, a 2 be taken 



as free parameters [ES 



Given the factorized Hamiltonian (|22|), one can then write down the amplitude for 
a decay. For example, if X~ is a meson made of valence quarks d and u, 

Amp(B° - D+X-) = ^V^V; fe a 1 (X-|M)L d |0)( J D + |(c6) hadM |5 ) (24) 
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where we have from Lorentz invariance 



(X~\(ud)^ d \0) = -ifxq" (for X: pseudo scalar) (25) 
(X~|(w<i)h ad |0) = fx m x^ (for X: vector or axial vector) (26) 

with fx being a parameter of energy dimension (called the decay constant). The current 
matrix element is the same as that appears in the corresponding semileptonic decay j|T| 
evaluated at q 2 = m 2 : 

(D+\(cb) had »\B°) = ( Pb + P D - ™*^q) F 1 (q 2 ) + ^^q,F (q 2 ) (27) 

where Fq and F\ are longitudinal and transverse form factors respectively [one can easily 
verify that the coefficient of F\ satisfies {...)^ = 0]. For the case of pion emission, the 
transverse component exactly vanishes (by definition) and we have 

Amp(B° - D+Tr-) = -i^V* d V cbai f n (m 2 B - m 2 D )F (ml). (28) 

The form factors F 0j i may be either obtained from semileptonic decays or calculated by 
models such as the relativistic harmonic oscillator model together with the pole dominance 
41fl . They are relatively slowly varying functions of order 1. In addition, the heavy quark 



effective theory allows us to relate all form factors for transitions between heavy mesons 



to a universal form factor ||42| . Similar procedures are applied to other decay modes. 

In general, we may distinguish three classes of decays when we consider two-body 
decays of heavy mesons mediated by operators of the types Oi j2 in spectator mode (i.e. 
the light quark in the parent meson does not participate in the weak decay) Q : 

Class 1 Only the first term in ( |2"2"D contributes and the amplitude is proportional to a\\ 
e.g. B° -> D+iT-. 



Class 2 Only the second term in (^) contributes and the amplitude is proportional to 
a 2 ; e.g. B — > _D°7r°. Sometimes called 'color-suppressed' decays. 



Class 3 Both terms in ( p22|) contribute and the amplitude contains both a\ and a 2 ; e.g. 
B~ -> D°7r- 

Some comments are in order. If both final-state particles are charged, then it is Class 
1, if both are neutral, then it is Class 2, if one is neutral and the other is charged, then 
it depends. In B — > D + n~, the current B — > D emits a n and thus the pion decay 
constant f n is involved. In B —>■ D°tt , the current B — > 7r emits a D meson and thus 
the D meson decay constant fr> is involved. In B~ — > D°tt~, a class 1 amplitude and 
a class 2 amplitude interfere and thus both f n and fjj are involved. Also, note that in 
B — > D°tt°, the 'color transparency' argument does not apply since the color-singlet pair 
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Figure 6: The annihilation and VT-exchange processes. 



passing through the cloud is now cu pair which are moving quite slowly, and it may form 
a D meson before leaving the cloud. Thus, factorization may not be a good assumption 
in this case. 

Heavy mesons may also decay through valence quark annihilation or ^-exchange 
processes as shown in Figure |6] which are also mediated by interactions of types Oi 2- 
Such processes have been discussed in the context of the lifetime difference between D + 
and D°, but thought to be helicity-suppressed jKJ, and also suppressed by form factor 
effect when two-body decays are considered ||46|[ . It was suggested, however, that the 
helicity suppression may be lifted when soft gluon effects are taken into account [ pT7fl . 
Even though annihilation/exchange processes are usually ignored in B decays, it has not 
been proven that they do not significantly contribute in all types of decays. 



3.3 Experimental Test of Factorization 

The decays B — ► PP, PV have definite final spin state, where P is a pseudo scalar meson 
and V a vector meson, thus the decay rate is the only dynamical parameter that can be 
tested. On the other hand, the decays B — > VV has three possible helicity amplitudes 
which can also be compared against prediction of factorization. 

For the test of decay rates, we take B — ► D* + X~ with X~ being 7r _ , p~ , or 
di . As described above, factorization allows us to estimate the decay rates of these 
modes from the g 2 -dependent form factors of the corresponding semileptonic decay B — > 
D* + l~v. In other words, there is a simple relation between the differential decay rate of 
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the semileptonic mode at q 2 = m 2 x and the corresponding non-leptonic decay rate, which 
can be conveniently written as pc 



Br(B° -> D* + X- 



^fx\VudY (29) 



where fx is the decay constant of the meson X. No QCD correction is included in the 
expression on the right hand side [fi9|| . If QCD correction is to be included, a reasonable 
choice would be to add (Ci + C 2 /3) 2 to the right hand side of (|2"9"|). This is because, in fl2~T|), 
the contribution from the octet current has been shown to be suppressed in the decays in 
question as shown by Dugan and Grinstein |35fl . However, C\ + C2/3 is unity to the first 



order due to the relation C\C\ = 1; thus, we will proceed without QCD correction. The 
above formula is applicable for X being any spin-1 particle or any light spin-0 particle 
(assuming factorization, of course) [f£8| . When the particle X is spin-0, it cannot replace 
all the helicity degrees of freedom of the D* appearing in the semileptonic decay, and the 
formula is correct only in the limit of nix <C ttib- The correction for pion, however, is 
negligible (~ 0.5%). If X is spin-1, then no such restriction applies. If D* is replaced by 
D, then a similar helicity projection factor should be included. 

The procedure of the test is to measure the decay rate B — > D*X and the differential 
semileptonic rate dT/dq 2 to obtain the ratio R, and then compare it to the value expected 
from factorization: Gtt 2 fx\Vub\ 2 ■ The q 2 distribution of the semileptonic decay B — > 
D* + l~v is shown in Figure [7], which is a combination of ARGUS J5(| and CLEO |5l| data. 
The shape is fit to three different models PH [52], Q to obtain the value at given q 2 . 



The decay constants can be obtained by the leptonic decay rate [54 




r(7r- - irvj = ^m^mjf 1 - ^ ) , (30) 
87r ^ \ m 2 

for pion which gives f n = 132 MeV. From the tau decay rates 

i 

167T r ' J " y m 2 J \ m' 2 
where V is a vector or axial vector, we get f p = 197 ± 3 MeV and f ai = 178 ± 28 MeV. 



r(r" - V-V T ) = ^-G 2 F \V ud \ 2 f 2 ( 1 - ^ U + (31) 



Including the effect of decay width of meson ||56fl , these go up to f p = 210 ± 3 MeV and 
f ai = 201 ± 32 MeV. Using the isospin symmetry relation f p - = f p o, the decay constant 
of p can also be obtained from T(p° —> e + e~) measured in e + e~ — > p° by 

r(V° - e + e-) = ^cyf 2 (32) 



where cy = 1/2 for p this gives / p = 216 ± 5 MeV which we will use. 

Table ^| summarizes the result of the comparison. Note that in taking the ratio (|29|) , 
uncertainty in D* detection efficiency is canceled. This of course assumes that same D* 
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Figure 7: The distribution of the lepton-neutrino invariant mass (q 2 ) for the process 
B -> D*+l~v as measured by CLEO and ARGUS. 

and D branching ratios are used for the measurements of D*lv mode and D*X mode; 
a correction has been made to the values of Figure [F] using the new measurements from 
CLEO pT], |53fl . The agreement is quite good for n and p. For ai, the measured R 
is about a factor of two larger than the expected value, but statistically it is only 1.5 
sigma's. This could well be due to breakdown of factorization at ai mass of 1.26 GeV. 
The branching ratio of D*ai is determined assuming that the D* + 7t + 7t~tt~ final state 
with 1.0 < Ms n < 1.6 GeV is dominated by a±. Figure § shows the 3tt invariant mass 
distribution for the decay mode B — > D* + tt + 7t~tt~ . The a\ peak is clearly seen, and 
amount of non a\ contribution is quite small. 

As stated earlier, for B — > VV decays there are helicity degrees of freedom which 
cannot be uniquely determined by kinematics. The factorization assumption leads to 
specific prediction for helicity amplitudes which can then be tested experimentally. For 
example, once the matrix element is factorized as 

Amp(B - D*p) oc (p\(ud)UU)(D*\(cb) h ^B), (33) 

then by Lorentz invariance the rho production part (p|(M^)h a d|0) is proportional to the p 
polarization vector e M [see (p6|)] . It then acts the same way as the polarization vector of W 
in semileptonic decay resulting in the same p polarization as that of the W in semileptonic 
decay at q 2 = m 2 . If factorization is not valid, this argument cannot hold, and thus it 
serves as a good check of factorization. 
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The polarization of p can be measured by the distribution of p~ — > 7r°7r~ polar decay 
angle 9 P in the p rest frame with respect to the D* direction in the same frame. Longitudi- 
nal polarization (helicity=0) would have cos 2 9 P distribution while transverse polarization 
(helicity=±l) would have sin 2 9 P distribution. Or equivalently, one can measure the decay 
angle of D* (9d) in the same way since the helicity of D* is the same as that of p. In fact, 
the angular distribution can be written as 

oc — — - sin 2 9 P sin 2 9b + cll cos 2 9 p cos 2 6d (34) 



d cos 9 p d cos Op 
where is the fraction of longitudinal polarization 





#0 


2 


\H + 


2 + 


H 


2 + 


H 


2 



a L ~ FT? ,2 I \tt 12 I IZT 12 ( 35 ) 



with being the three helicity amplitudes. Figure |^ shows the distributions of 9 P 

and 9b for data. A simultaneous fit to the two angles gives 

a L = 0.93 ±0.05 ±0.05 (CLEO) (36) 

where the first error is statistical and the second systematic which includes uncertainty in 
background subtraction and detection efficiencies. The experimental measurement of the 
polarization in the semileptonic decay is unfortunately not available at this point, and we 
have to compare the above measurement to absolute theoretical prediction which requires 
some assumption on form factors. One estimate using HQET jP| gives 

a L = 0.88 (factorization + HQET) (37) 

which is in agreement with the data. 



X 


Br(D* + X-) dBr/dq 2 ^(measured) 
(%Y a ^ (%/GeV 2 ) (GeV 2 ) 


fx R = 6n 2 \V ud 
(GeV) (GeV 2 ) 


2 £2 

Jx 


TT 
P 


0.265±0.036 0.23±0.05 1.15±0.30 
0.735±0.106 0.25±0.04 2.94±0.63 
1.32 ±0.30^ 0.32±0.04 4.13±1.07 


0.132±0.0005 0.98±0.01 
0.216±0.005 2.63±0.12 
0.201±0.032 2.27±0.72 



(a) The errors are statistical only. 

(b) It is assumed that D* + a^ dominates D* + 7t + 7t~tt~ mode where the 3n mass is between 
1.0 and 1.6 GeV. 



Table 3: Test of factorization. Branching fraction of B — > D*X is compared to the 
corresponding semileptonic decay evaluated at q 2 = m x . 
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Figure 8: The 3n mass distribution in the decay B° — > D* + n + n n . (a) Monte Carlo 
simulation for B° — > D* + di,a{ — > p°7r _ . (a) Monte Carlo simulation where p°7r~ is 
uniform in phase space, (c) Data with S-mass side bands subtracted. 



Figure 9: The angular distributions for D* decay angle (a) and the p decay angle (b) in 
B° -> D* + p~. 
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This helicity=0 dominance can be intuitively understood as follows (see the figure 
above). When the ud pair is emitted, they are nearly collinear, and the helicities are 
left-handed for the d quark and right-handed for the u quark. Therefore the total helicity 
is zero which is transferred to the final p meson assuming that there is no final state 
interaction that changes the spin state. This feature is independent of specific choice of 
form factors, while it does assume factorization. 



3.4 Extraction of a\ and a% 

In this section, we will take the coefficients a% and a 2 as free parameters in the factorized 
effective Hamiltonian (^3), and try to find their values by fitting to measured branching 
ratios. First, we will use B — ► Dti and ipK decays to demonstrate the procedure, then a 
global fit to clean modes will be performed. 

In order to extract a\ and a-2, we need the form factors of B — ► D transition. This 
is quite well known; we will use the result of the fit to the universal form factor under 
the framework of HQET. The relevant value here is Fo(q 2 = m 2 ) = 0.58. For the B — > 7r 
or K transition, there is no experimental data, and a model calculation is used where 
the overlap of B and the light meson wave functions is obtained by relativistic harmonic 
oscillator model and the q 2 dependence is given by pole dominance. The coefficients of 
012 below are taken from Reference |27 . 



Class-I (determination of ai): The decay amplitude of B° — ► D + n (or for any 



two-body decay B — > PP) is given by 



8nM 



'' \Amp\ 2 (38) 



2 



where p is the momentum in the B rest frame. Using the factorized amplitude ( |28| ) 
together with V cb = 0.045, V ud = 0.975, G F = 1.166 x 10" 5 (GeV~ 2 ), F (ml) = 0.58, and 
r B = 1.18 ps, we get 

Br(D + n-) = 0.264a 2 (%). (39) 

The measured branching ratio is Br(D + iT~) = 0.29 ± 0.04% from CLEO, where the error 
is statistical only. It then gives a\ = 1.1. 

Class-II (determination of a 2 ): In B° — > D°ti , D meson is emitted and the transition 
is from B to 71. Proceeding the same way as before, we get 

5r(DV) = 0.20l(M^) 2 a 2 (%). (40) 

where the isospin factor 1/2 is included (n° is half uu and half dd). Experimentally, 
only upper limit exists for this mode: a recent number from CLEO is Br(D°TT°) < 
0.035%(90%C.L.), which corresponds to \a 2 \ < 0.4. 

The decay B~ — > ipK~ is also a Cabbibo-favored Class-II decay. The Mb distribu- 
tion by CLEO is shown in Figure |10|(a) together with other related modes that are used 
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Figure 10: The beam-constrained masses for B 
ipK*-(c), and B° -> ipK*°(d) (CLEO collaboration) 



in the combined fit later; there are about 60 signal events with little background. Mass 
peaks for B~ — > ipK~ and ipK*° by CDF are shown in Figure 11 [59|. The transition is 
B —>■ K and if> is emitted. The decay constant of if) can be obtained from its e + e~ width: 
= 384 ± 14MeV. The expected branching ratio is 



Br{ipK-) = 1.819e^ (%) 



(41) 



where the large coefficient is primarily due to the large decay constant of if). The measure- 
ment Br(ipK~) = 0.110 ± 0.015 (CLEO) gives \a 2 \ = 0.26. One point of caution is that 
ci2 in b — > ccs transition is likely to be different from a 2 in b — > cud transition. In fact, the 
values of themselves are expected to be different as seen in (|lg|). Nonetheless, they 
are often assumed to be the same and we will proceed with this assumption for now. 

Class-Ill (determination of g^/g^): As stated earlier, for Class-II and Class-Ill de- 
cays, the factorization assumption is not well founded. However, if we assume the fac- 
torized Hamiltonian (|22|), we can obtain the sign as well as the absolute value of 02/01 
through the interference of the two types of diagrams shown in Figure f|. For example, 
the branching fraction of B — > D°tt~ (normalized to B — > D + tt~) is given by 



Br{D 71 
Br(D+n-) 



aA 220 



1 2 



(42) 
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Figure 11: Invariant mass peaks for B — > ipK (a) and B° — > ipK*°(b) by the CDF 
collaboration. 

The ratio measured by CLEO is 1.84 ±0.24 ±0.29, and this leads to a 2 /a 1 = 0.29 ±0.11. 
The positive sign is a direct consequence of Br(D°7T~) > Br(D + 7i~). 

Tables |1B] summarize measurements and expected branching ratios from the factor- 



ization model as calculated in Reference [27|. The agreements are excellent in all cases 



In order to obtain more accurate value for a\ we fit four Class-I modes, B — > 
D + tt~, D + p~, D* + n~, and D* + p~ . For a 2 , we use the Class-II modes B° — > ipK ,if)K* 
and B~ — > ipK~ ,ipK*~ . The result is 

|ai| = 1.15 ±0.04 ±0.04 ±0.09, \a 2 \ = 0.26 ± 0.01 ± 0.01 ± 0.02 (43) 

where the first error is statistical, the second and the third are systematic. The third error 
is due to the uncertainty in the ratio of production and that of lifetimes of charged vs 
neutral B mesons. The relevant quantity is (/+t+)/(/o t o) where / + , f are the production 
fractions and r + , tq are the lifetimes. This value is sometimes assumed to be unity. A 
measurement from Br(B~ — > D*°lu) / Br(B° — > D* + lu) |6^] is 



= 1.2 ±0.20 ±0.10 ±0.16. (CLEO). (44) 

Jo ?o 

For determination of 02/01, we use the following four ratios of branching fractions: 
B(D°it-)/B(D + it-), B(D°p-)/B(D + p-), B(D*°n-)/ B(D* + ir~), and B(D*°p-)/B(D* + p- 
to obtain 

— = 0.23 ±0.04 ±0.03 ±0.10 (45) 
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D 


CLEO (%) 


ARGUS (%) 


Model (%) H 


ai = 1.15 


D + 7T~ 

1 J l\ 


0.29 ± 0.04 ± 0.03 ±0.05 a 


0.48 ± 0.11 ±0.11 d 


0.264a? 


0.35 


D+p- 


0.81 ±0.11 ±0.12 ±0.13 a 


0.9 ± 0.5 ± 0.3 d 


0.621a? 


0.82 


D*+TT- 


0.26 ± 0.03 ± 0.03 ±0.01 a 


0.28 ± 0.09 ±0.06 d 


0.254a? 


0.34 


D* + p- 


0.74 ± 0.10 ± 0.13 ±0.03 a 


0.7 ± 0.3 ± 0.3 d 


0.702a? 


0.93 


D* + a^ e 


1.26 ± 0.20 ± 0.14 ±0.04 a 




0.97a?(/ al /0.22) 2 


1.28 


n**+ 

-^(2460) n 


< 0.18 a 








D + D~ 


1.2±0.7 fe 


1.7 ± 1.3±0.6 C 


1.213a?(/ Ds /0.28) 2 
0.859a?(/ Ds ,/0.28) 2 
0.824a?(W0.28) 2 
2.203a 2 (/ Ds */0.28) 2 


1.60 


D + D*- 




2.7 ± 1.7±0.9 C 


1.14 


D* + D- 


2.4 ± 1.4 fe 


1.4 ± 1.0±0.3 C 


1.09 


£>*+£>*- 




2.6 ± 1.4±0.6 C 


2.91 










2.9 ± 1.3 b 


2.4 ± 1.2±0.4 C 


1.215a?(W0.28) 2 
0.862a?(/zW0.28) 2 
0.828a?(/ Ds /0.28) 2 


1.61 


D°D*~ 




1.6 ± 1.2 ± 0.3 C 


1.14 


D*°D; 




1.3±0.9±0.2 C 


1.10 


D*°D*- 




3.1±1.6±0.5 C 


2.206a?(W0.28) 2 


2.92 



)7T 



a. Preliminary result to be submitted to Phys. Rev. D. The first error is statistical, the 
second systematic, and the third error is due to uncertainties of D branching ratios. 

b. Reference |5tJ, Br(Dg 

c. Reference [61], Br(Dg 

d. Reference [[J. 

e. All events with 37r mass between 1.0 and 1.6 GeV (after background subtraction) are 
assumed to be a\. 



07T 



2% is used. 
2.7% is used. 



Table 4: Class-I Branching ratios 
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CLEO (%) 


ARGUS (%) 


Model (%) [0 


ax = 0.26 


D°n° 


< 0.035 a 




0.201a2(/ D /0.22) 2 
0.136a 2 (/ D /0.22) 2 
0.213a^(/ £) */0.22) i 

\J.ZZ0CL2\J D*l U.ZZJ 


0.014 


D°p° 


< 0.042 a 




0.009 


D* U 7T U 


< 0.072 a 




0.014 


n*o „o 
u p 


<s u.uyz 




n m k 

U.UIO 


D°r] 


< 0.075 a 








D°7]' 


< 0.074 a 








D°LU 


< 0.048 a 








D*°T] 


< 0.086 a 








D*°ri' 


< 0.36 a 








D*°LU 


< 0.13 a 








iU iV 


0.075 ± 0.024 ±0.008 a 


0.08 ± 0.06 ± 0.02 b 


1.817a 2 . 


0.123 


iblC* 

iU XV 


0.169 ± 0.031 ±0.018 a 


0.11 ± 0.05 ±0.02 b 


2.927a 2 


0.198 




< 0.08 a 


< 0.28 fe 


1.065a 2 


0.072 




< 0.19 a 


< 0.23 fe 


1.965af. 


0.133 




< 0.27 a 










< 0.21 a 
















0.110 ± 0.015 ±0.009 a 


0.07 ±0.03 ±0.01" 


1.819a 2 


0.123 


ipK*~ 


0.178 ± 0.051 ±0.023 a 


0.16 ± 0.11 ±0.03 6 


2.932a 2 


0.198 


fK- 


0.061 ±0.023 ±0.015 a 


0.18 ± 0.08 ±0.04 b 


1.068a 2 


0.072 


ijj'K*~ 


< 0.30 a 


< 0.49 6 


1.971af. 


0.133 


XdK- 


0.097 ± 0.040 ±0.009 a 








XciK*- 


< 0.21 a 









a. Preliminary result to be submitted to Phys. Rev. D. 

b. Reference ||. Modes involving a K s are multiplied by two to obtain the branching 
ratios for K° . 



Table 5: Class-II Branching ratios 
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B~ 


CLEO (%) 


ARGUS (%) 


Model (%) |7J 


CL\ = 1.15 

a 2 = 0.26 


D°n- 


0.55 ±0.04 
±0.03±0.02 a 


0.20 ± 0.08 ±0.06 b 


0.265(ai + 1.230a 2 (/ D /0.22)) 2 


0.57 


D°p- 


1.35 ±0.12 
±0.12±0.04 a 


1.3±0.4±0.4 6 


0.622(ai + 0.662a 2 (/ D /0.22)) 2 


1.09 


D*°tt- 


0.49 ±0.07 
±0.06±0.03 a 


0.40 ± 0.14 ±0.12 b 


0.255(ax + 1.292a 2 (/W0.22)) 2 


0.56 


D*°p- 


1.68 ±0.21 
±0.22±0.08 a 


1.0±0.6±0.4 6 


0.703[a 2 + 0.635a 2 (/ D ,/0.22) 2 
+1.487o 1 a 2 (/u*/0.22)] 


1.27 


D*°a^ c 


1.88 ±0.40 
±0.30±0.10 a 








-^(2420) 71 


0.11 ±0.05 
±0.04 ±0.03" 








-^(2460) 71 


< 0.15 a 








D (2420) P 


< 0.14 a 








D*U60)P 


< 0.5 a 









a. Preliminary result to be submitted to Phys. Rev. D. The first error is statistical, the 
second systematic, and the third error is due to uncertainties of D branching ratios. 

b. Reference [§]. 

c. All events with 37r mass between 1.0 and 1.6 GeV (after background subtraction) are 
assumed to be a\. 



Table 6: Class-Ill Branching ratios 



23 



where (/+T + )/(/o r o) = 1-2 is used and the last error is due to the uncertainty in this 
quantity. The absolute value of 02/01 is consistent with the value obtained above which 
is 0.26/1.15 = 0.23, and the negative sign seems to be excluded. From fl2~3f) , we have 



a x ~ C 2 + £d ~* * ~ 1 - (C a /Ci)(o a / ai ) ' 1 j 

Using Ci = 1.11, C 2 = —0.26, the negative value a 2 /oi = —0.23 corresponds to £ = 0.01 
and the positive value 02/01 = 0.23 corresponds to £ = 0.44. Thus, £ = as suggested by 
an analysis of charm decays |28| seems to be excluded in the B decays. However, one has 



to keep in mind that in the analysis above, the factorization was applied to questionable 
cases where emitted meson is heavy. Also, the factorization is not expected to hold well 
for charm decays, so the formulation using 01^ itself is in question in charm decays. 

So far only Class-II modes observed are for b — *> ccs only. As one can see from the 
table, however, the present sensitivity is close to the expected values for the D°7r° and 
related modes. It is likely that these modes will be observed soon. 

3.5 Final State Interaction 

The factorization assumes that effect of final state interaction is negligible. Therefore any 
test that is sensitive to final state interaction is also a test of factorization. 

One way is to perform an isospin analysis on a set of isospin-related modes. For 
example, the relevant Hamiltonian for B — > Dtx decays has isospin structure \I,I Z ) = 
|1, — 1) (i.e. b — ► cud - simply a creation of ud pair as long as isospin is concerned). 
Separating the Hamiltonian to an isospin violating part S and an isospin conserving part 
h, and treating S as if it is a particle (spurion), we can write, for example 

SB° = y||3/2,-l/2)-y||l/2,-l/2) 

D + n- = ^|3/2,-l/2)-^||l/2,-l/2). (47) 

Applying similar isospin decomposition to SB~, B° — > D°7i°, and B~ — > D°tt~, we get 

Amp(SB° -> D+U-) = i(3/2,-l/2|o|3/2,-l/2) + ^(l/2,-l/2|o|l/2,-l/2) 

Amp(SB° D°7T°) = — (3/2,-l/2|o|3/2,-l/2) - — (1/2, -l/2|o|l/2, -1/2) 

3 3 

Amp(SB- -> D°7T) = (3/2, -3/2|/i|3/2, -3/2) (48) 

Because of isospin invariance of h, the matrix elements depend only on the magnitude 
of the isospin: (3/2, -3/2|o|3/2, -3/2) = (3/2, -l/2|/i|3/2, -1/2) which we define to be 
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y/SAa. Together with a definition (1/2, -1/2 1 to 1 1/2, -1/2) = we then obtain 

/T /o" 

Amp(D + 7T~) = 




Amp(D°n°) = J^As-^Ai (1.9) 
v4mp(D 7T~) = 73^3 



where A3 and Ai are the isospin 3/2 and 1/2 amplitudes respectively. There are three 
unknown parameters: | ^4.3 1 , \A±\, and 5 = a,rg(A3/Ai). Since there are three measure- 
ments of decay rates, one can solve for the three unknowns. Then the non-zero phase 
8 signifies the existence of final state interaction. Unfortunately, the D tt mode is not 
observed yet at this point; we expect, however, that it will be observed sometime soon as 
mentioned earlier. 

One could go further along this line if one is bold enough. One can set 5 = 
and recalculate the decay rates that would have been without the final state interaction. 
Then those rates may be compared with what is expected by factorization. In fact, a 
phenomenologically successful analysis of charm decay was performed in such manner 
f28|. However, there is no guarantee that all the effect of final state interaction can be 
taken away by this method. There may be interactions with other final states, for example. 

Another possibility is to look at the azimuthal angular distribution in B — > VV 
decays. Taking B — > D*p as an example, the angular distribution is given by 

cc (\H + \ 2 + \H_\ 2 )s 2 D s 2 p + A\Ho\ 2 c 2 D c 2 p 

+2Rq(H* + H J) s 2 D s 2 p cos 2 X + 2\m(H* + H J) s 2 D s 2 p sin2 X (50) 
+4Re(if*if + Hl_H ) s D c D s p c p cos x + 4Im.(H+H - H*LH )s D c D s p c p sinx 

where 9d, p are the polar decay angle of D and p decays as before, and x is the azimuthal 
angle between the two decay planes. We have used a short hand: cp = cos#d, sd = sin^ 
etc. If there is no final-state interaction and there is no CP violation, then all the helicity 
amplitudes are relatively real. The effect of CP violation would show up as difference of 
angular distribution (as well as difference in total decay rate) between B and B decays 
|k|. For Cabbibo-favored modes such as D*p, we do not expect significant CP violation. 
Thus, existence of terms proportional to sin x or sin 2% signals final state interaction |M 



This analysis should be able to be done with dataset presently available, but thus far not 
completed. 



4 Suppressed Decays 

Now we move to rare decays which are typically Cabbibo-suppressed. We start from 
charm-less two-body decays. 
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4.1 B Decays to Two Charmless Mesons 

Each of the processes ~7i + , 7r + 7r could proceed through two types of diagrams: 

spectator and penguin (Figure 0). When there exist more than one diagram with dif- 
ferent weak interaction phases and different final state interaction phases (i.e. strong 
interaction), there can be CP violating decay asymmetries [65[] as seen below. Suppose 
two diagrams contribute to a decay B — > / with amplitudes A\ and A<i& & where A% t 2 
are the weak amplitudes and 5 is the FSI phase difference. Since only relative phases 
matter, the weak and strong phases of the first diagram are assumed to be zero. For the 
corresponding B — > f decay, the weak phase changes its sign but the strong phase does 
not. This leads to a decay asymmetry: 

Amp(B -+f) = A x + A 2 e iS , Amp(B -+f)=A{ + A* 2 e iS (Ax : real) (51) 

(52) 



In our case, the weak phase of each diagram is given by that of the CKM matrix elements 
which multiply the entire amplitude as coefficients. Thus we expect that there is a weak 
phase difference as shown in the figure. The strong phases, however, are difficult to 
estimate. 

If we assume the flavor SU(3) symmetry, then the ratio of amplitudes are 



Kit 



TXTX 



Ktt 

A 



~ T (53) 

penguin A 



spectator 

where A is the Cabbibo suppression factor (~ 0.2). It is expected that the spectator 
diagram will dominate in B° — > 7r + 7r _ . Then if there is no penguin contribution, the 
K~n + branching ratio should be A 2 ~ 0.04 times smaller than that of 7r + 7r~. Thus, if the 
rate of K~tt + is comparable or greater than 7t + tt~, then it is likely that the K~tt + rate 
is dominated by the penguin diagram. When there is a large disparity in magnitudes of 
the two diagrams, the expected CP violation will be small independent of the phases. 

One should note, however, that there is a possibility that B — > Ktc can occur 
through final state re-scatterings. This could occur through intermediate states involving 
two charmed mesons as 

B° -> D~Dg -> K+tT, B° -> D-D + -> it-it + (54) 
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Figure 12: Diagrams that can contribute to B — > Kir, 
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which corresponds to replacing the top quark loop in the penguin diagrams by a charm 
quark which will be on-shell as shown below and can be considered to be a dispersive 
version of penguin diagram. 

(55) 



Such process will result in a large FSI phase, and can interfere with the top quark penguin 
diagram to generate a CP violation as originally postulated by Bander, Silverman and 
Soni . 

Approximate rate of tt + tt~ can be estimated from the measured B° — > D~tt + rate 
quite reliably: 

V 2 

"" Br{D + 7i-) ~ 1 x 10" 5 (56) 



Br (71 7T 



V cb i 

where the effect of form factor will reduce it somewhat and that of phase space will 
increase it somewhat. The estimation of the Kit rate requires the coefficient of the 
penguin operator, and the uncertainty is greater; the theoretical estimates are in the 
same range as the tttt mode |)7] . 

Experimentally, the signature on T(4S) is a rather spectacular high- momentum 
back-to-back tracks of p ~ 2.6 GeV. This is the maximum momentum a -B-decay can 
emit and the background is dominated by continuum events; thus, cuts are made on 
event shapes to reject 2-jet like events and the fast back-to-back tracks are required not 
to be aligned with the jet axis of event. For a BB pair event, the event shape is spherical 
and there is little correlation between the event axis and the direction of the back-to-back 
tracks. Then, as before, the energy difference AE and the beam-constrained mass Mb is 
used to select the candidates [see (|3[)]. 

When masses are correctly assigned to the tracks, the AE resolution is 25 MeV. 
The dE/dx information in the drift chamber is used to separate kaon and pion. The 
dE/dx resolution is 6.5% and provides 1.8a K — tt separation per track. Each candidate 
is assigned the most likely masses (tttt, Ktt, or KK), then AE is calculated. The beam- 
constrained mass, on the other hand, does not depend on the mass assignments and the 
resolution is 2.5 MeV. Figure |l3](a) shows the Mb distribution for Ktt and tttt candidates 
after 2-er cut on AE around zero. The shaded events are the tttt candidates. One can see 
an enhancement at the nominal B mass of 5.280 GeV. The AE distribution after the 2-er 
cut on Mb is shown in Figure |T3|(b). Again, there is a peak around the nominal region 
near AE = 0. For the final extraction of numbers, an un-binned maximum likelihood fit 
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Figure 13: Sum of Kti sample and tttt sample, (a) The beam-constrained mass Mb after 
the 2-cr cut on AE. (b) AE distribution after the 2-a cut on Mb- The shaded events are 
the events assigned to be tttt. 

is performed with AE, Mb, dE/dx, and an event shape variable as parameters. Here AE 
is calculated assuming tttt. The result is shown in Table [7| |68[] . 

When AE is calculated assuming tttt, the value shifts down by 42 MeV if the actual 
tracks are Ktt. Since the AE resolution is 25 MeV, this by itself can provide 1.7a sep- 
aration between Ktt and tttt. The available particle identifications are not good enough 
to cleanly separate the two. When Kti and tttt are combined there is a substantial signal 
of about 3.5a. The central value of 7r + 7r~ mode is consistent with the expected value of 
1 x 10 5 . If we take the central value of the K + tt~ mode at its face value, then penguin 
diagrams (t-loop or the re-scattering c-loop) are likely to be dominating the K + n~ mode. 



Mode 


Br(l(r 5 ) 


Upper Limit (10 5 ) 




I.SHq'6 ± 0.2 


2.9 


K+ir- 


1.1±°; 6 7 ±0.2 


2.6 


K + K- 




0.7 


K + 7T~ + 7T + 7T~ 


2.41^ ± 0.2 



Table 7: Measured branching fractions and 90% confidence level upper limits. 
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4.2 b to s Radiative Decays 



Another rare process a penguin diagram is expected to contribute is the radiative b 
transition through emission and re-absorption of W. 



(57) 



At the lowest order, the GIM suppression is operative and it depends on the top mass 
(nit) strongly, and BriB —>■ X s / y) changes from 0.5 x 10~ 4 at m t = 100 GeV to 1.4 x 10~ 4 
at m t = 200 GeV. With QCD correction fl6H|] , the GIM suppression is loosened ('soft' GIM 
suppression) and as a result the rate is substantially enhanced and becomes a slow function 
of mt- The enhancement factor is ~5 at m t = 120 GeV to give BriB — > X s ^y) ~ 3.5 x 10~ 4 . 
Theoretical estimate for the exclusive mode B — ■> K*^y is more uncertain due to the 
unknown transition matrix element B —>■ K* |70|| . One estimate based on HQET gives 



Br{B -»■ K*j) = (1.4 - 4.9) x 10~ 5 [[TT 



The experimental signature [7^] is a monochromatic hard photon (2.6 GeV) recoiling 
against K* -> Kir decay. We look for both B° -> K*°-f and B~ -> K*-^. The iT*'s are 
searched for in the modes K*° — > K + tt~ and K*~ — > K~tt°, K s tt~. Again the background 
is dominated by continuum events since such high-energy photon is at the kinematic limit 
of B decay. The continuum backgrounds are reduced by requiring that the events be not 
2-jet like and that the hard photon be not aligned to the event axis. If the photon forms 



a 7T° or rj with another photon then it is rejected. Figure [TJ] shows the Mb distribution 



after the cut \AE\ < 90 MeV (2.2cr). There is a clear signal observed with 6.6 ± 2.8 
events in B° mode and 4.1 ± 2.3 events in B~ modes. The branching fractions are 

Br(B° — > K*°^y) = (4.0 ± 1.7 ±0.8) x 10~ 5 , 

Br(B~ ->• K*~i) = (5.7 ±3.1 ± 1.1) x 10~ 5 (CLEO) (58) 

If we assume isospin symmetry, then 

Br(B° -> K*° 7 ) = Br(B~ -> K*--f) = (4.5 ± 1.5 ± 0.9) x 10" 5 (59) 

which is consistent with theoretical estimates based on the standard model where the 
penguin contribution dominates. Another possibility is that the B — > K*^ transition 
may occur through ipK* by vector dominance |73[ 



B — > ipK* — > '-/K* (vector dominance). (60) 
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Figure 14: The Mb distribution for B 



if* 7 after AE cut. 



Figure 15: Single photon spectrum after continuum subtraction. The signal b — * 57 would 
show up in the region 2.2 to 2.7 GeV. 
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or other long distance effects WM. Such processes have been estimated and found to be 



at least an order of magnitude smaller than the observed rate. 

The inclusive transition B — > X s ^ can be searched by looking for the hard photon 
without reconstructing X s where the mass of X s lies in the typical strange meson region 
(0.5 to 2 GeV). Similar cuts as before to reduce continuum backgrounds are applied. 
Figure [15| shows the continuum-subtracted (see Section 1) photon spectrum. The signal 
region is around 2.2 to 2.7 GeV. There seems to be some enhancement, but it is not 
statistically significant; thus, we set an upper limit 

Br(b -> sj) < 5.4 x 10" 4 (CLEO[75]). (61) 

Such measurement places stringent constraints on non-standard physics, in particular two- 
Higgs-doublet models |76[| . The W-top loop can be replaced by loops involving charged 



Higgs, neutralinos, gluinos, and squarks etc |77|. For example, in the minimal supersym- 



metric model with two Higgs doublets, the mass of the CP-odd neutral Higgs A (which 
is related to the mass of the charged Higgs) is ruled out for m^o < 250 GeV. ff8| . 
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